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ABSTRACT 

Trade-off  models  are  developed  for  siting  inventories 
of  repair  parts  by  Naval  Supply  Center  (NSC)  which  must  sup¬ 
port  a  local  Naval  Air  Rework  Facility  (NARF) .  Two  strategies 
are  considered;  siting  at  the  NSC  with  direct  delivery  to  the 
NARF  and  siting  at  the  NARF.  Three  direct  delivery  alterna¬ 
tives  which  include  both  scheduled  and  unscheduled  delivery 
schemes  are  modeled  when  siting  is  at  the  NSC.  The  measure 
of  effectiveness  for  all  alternatives  is  the  expected  total 
costs  per  time  period.  Cost  elements  include  delivery  costs 
and  production  delay  costs.  Algorithms  for  solving  the 
trade-off  models  are  also  presented. 
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I.  INTRODUCTION 


In  the  consolidation  of  wholesale  supply  support 
between  the  Naval  Supply  Centers  (NSC)  at  Oakland,  San  Diego, 
and  Norfolk,  and  the  neighboring  Naval  Air  Station  (NAS) 
supply  centers,  the  question  of  supply  support  of  the  local 
Naval  Air  Rework  Facilities  (NARFs)  is  of  major  concern.  [1]  A 
goal  of  these  consolidations  is  not  to  degrade  service  to  the 
NARF  and  hence  the  question  of  whether  special  inventories 
should  be  located  at  the  customer's  site  is  raised.  One  obvi¬ 
ous  advantage  would  be  that  the  response  time  of  the  supply 
system  would  be  quicker  because  the  travel  time  would  be 
shorter  from  the  on-site  system  than  it  would  be  if  the 
material  were  held  at  the  NSC.  The  transportation  costs  and 
customer  delay  costs  would  also  be  less.  The  disadvantage  is 
that  additional  costs  are  incurred  in  maintaining  these 
on-site  systems. 

The  answer  to  this  question  of  supply  support  lies  in 
the  results  of  a  trade-off  analysis  which  seeks  to  find  a 
balance  between  the  three  cost  components  mentioned  above. 

Such  a  balance  may  result  in  all  the  inventory  at  the  customer 
site  and  none  at  the  NSC,  or  all  at  the  NSC  and  none  on-site, 
or  some  at  each  location.  The  purpose  of  this  report  is  to 
develop  the  models  needed  to  conduct  the  trade-off  analysis 
for  the  all-or-none  situations  since  they  correspond  to  the 
current  philosophy  of  the  NSCs  towards  NARF  support . 


Figure  1  suggests  the  details  which  should  be  con¬ 
sidered  in  the  trade-off  model.  The  element  denoted  as  OSIS 
is  the  on-site  inventory  system. 


Branch  A  of  the  figure  represents  the  alternative  of  direct 
delivery  from  the  NSC  to  the  NARF.  Branches  B  and  C  represent 
the  alternative  of  a  periodic  replenishment  to  the  OSIS  from 
the  NSC  and  direct  delivery  as  each  demand  occurs  from  the 
OSIS  to  the  NARF,  respectively.  Branch  A  deliveries  may  not 
be  immediate  upon  receipt  of  a  demand  in  contrast  to  Branch  C. 
In  fact,  the  most  economical  approach  may  be  to  delay  delivery 
until  several  demands  have  been  received. 

Deterministic  Demand — To  develop  a  basis  for  understanding 
the  more  complex  probabilistic  direct  delivery  models  to  be 
presented  later,  we  consider  the  case  of  deterministic  or 
known  demand.  Suppose  that  a  truck  has  a  capacity  of  n 
units  of  an  item.  Suppose  also  that  C^  is  the  cost  of  a 
trip  by  the  truck  from  the  NSC  to  the  NARF.  If  the  truck 
makes  a  trip  as  soon  as  the  NSC  receives  and  processes  a  unit 
demand,  then  the  cost  of  shipping  the  unit  is  C_.  If,  how¬ 
ever,  the  truck  waits  for  k  units,  then  the  shipping  cost 
per  unit  is 


and  obviously  the  cheapest  unit  cost  occurs  when  the  truck  is 
filled.  That  is, 

CT 
n  * 

However,  while  the  truck  is  waiting  for  a  full  load, 
the  units  waiting  will  cause  delays  in  production  and  these 
delays  will  result  in  extra  costs  to  the  NARF.  To  model  these 
delays,  let  us  suppose  that  a  unit  is  needed  every  t  units 
of  time  because  of  the  repair  schedule  and  that  the  cost  of  a 
delay  over  t  for  one  unit  is  C^.  If  we  wait  for  k  units 
to  be  accumulated,  the  total  delay  cost  will  be 

Cn  k  [k-1] 

-2— 2 -  (2) 

To  confirm  formula  (2),  consider  Figure  2. 


If  we  wait  until  k  units  have  accumulated,  the  truck 


will  not  leave  until  time  [k-l]t.  The  first  unit  demanded 
occurs  at  time  =  0  and  is  delayed  [k-1]  periods  of  length  t, 
the  second  unit  demanded  occurs  at  time  t,  and  is  delayed 
[k-2]  periods,  the  third  unit  is  demanded  at  time  2t,  and 
waits  [k-3]  periods,  and  so  on.  Only  the  kt^1  unit  doesn't 
wait.  The  total  waiting  time  in  periods  is  therefore 

[k-1]  +  [k-2]  +  [k-3]  +  -  +  1  +  0  .  (3) 

The  sum  given  by  (3)  can  be  written  is  a  short  form  as 

,  (4) 


and  when  we  multiply  (4)  by  the  cost  per  unit  delayed 

one  period  we  get  formula  (2) . 

The  average  delay  cost  per  unit  is  then  obtained  by 
dividing  (2)  by  k;  that  is, 


cD[k-i] 

2 


(5) 


The  total  average  cost  of  shipping  and  delay  can  now 
be  written  as  the  sum  of  formulas  (1)  and  (5) : 


C  (k) 


c  c  [k-i] 

"F  +  ~2 - 


(6) 


Equation  (6)  can  also  be  viewed  as  the  average  cost  per  period 
because  demand  is  a  known  function  of  time. 


Random  Demand — Equation  (6)  was  derived  based  on  the  assumption 
that  one  unit  of  a  repair  part  is  demanded  every  t  time 
units.  This  demand  pattern  would  correspond  to  a  repair  part 
which  is  replaced  in  every  component  undergoing  overhaul.  If, 
however,  a  repair  part  is  needed  only  p  percent  of  the  time 
in  such  a  component,  where  p  <  100%,  then  such  a  part  may  not 
be  demanded  every  t  unit  of  time  and  we  must  therefore  con¬ 
sider  minimizing  the  expected  total  costs  per  period.  And 
since  it  may  not  be  desirable  for  the  truck  to  always  be  filled 
when  it  makes  a  trip,  we  need  to  consider  some  alternative 
delivery  strategies.  Those  which  seem  most  appropriate  are: 

1.  The  truck  makes  a  delivery  at  the  end  of  every  N 
periods  if  there  is  at  least  one  demand  during  the 
N  periods. 

2.  The  truck  makes  a  delivery  as  soon  as  demands  have 
accumulated  to  a  specified  number  K  not  exceeding 
its  capacity  n. 

s  t 

3.  The  truck  makes  a  delivery  in  the  (N-l)  period 
following  the  first  demand  received  after  the  last 
delivery. 

Alternative  1  corresponds  to  scheduled  deliveries. 
Alternative  3  is  a  variant  of  scheduled  deliveries  where 
time  =  0  corresponds  to  the  time  of  the  first  demand.  Alter¬ 
native  2  corresponds  to  unscheduled  deliveries. 

To  compare  these  alternatives,  we  need  to  develop 
expressions  for  the  expected  average  total  costs  per  period 

5 


i 


because  the  time  between  deliveries  may  be  different  for  each 
alternative  and  thus  the  number  delivered  over  a  given  time 
interval  can  be  expected  to  be  different.  Chapters  II,  III, 
and  IV  present  the  derivations  needed. 

Another  set  of  models  are  needed  to  complete  the 
picture  depicted  in  Figure  1.  These  are  the  models  which  form 
the  basis  for  comparing  the  expected  total  costs  per  period 
of  operating  an  on-site  inventory  system  with  those  costs 
derived  in  Chapters  II,  III  and  IV.  The  derivations  are 
presented  in  Chapter  V. 

Chapter  VI  combines  the  models  from  Chapters  II 
through  V  into  a  structure  such  that  the  trade-off  analyses 
can  be  made.  Included  in  that  structure  is  the  impact  of 
constraints  imposed  by  time  standards  and  truck  capacity. 

These  constraints  are  derived  in  Chapters  II,  III  and  IV  for 
the  direct  delivery  alternatives.  The  final  chapter,  Chapter 
VII,  summarizes  the  modeling  efforts  and  the  results  of  the 
analysis.  It  also  makes  suggestions  for  model  refinements. 
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Introduction — Alternative  1  assumes  that  a  scheduled  delivery 
will  be  made  of  all  the  units  of  an  item  needed  to  meet  the 
demands  which  have  accumulated  by  the  end  of  N  periods.  If 
none  have  accumulated  the  scheduled  delivery  will  be  cancelled, 
Probabilities — We  assume  that  only  one  unit  of  a  given  repair 
part  is  needed  by  a  component  undergoing  repair.  We  also 
assume  that  the  probability  that  it  will  need  to  be  replaced 
is  a  constant  denoted  by  p.  Therefore,  the  probability  dis¬ 
tribution  of  the  demand  for  x  units  cf  the  repair  part  dur¬ 
ing  an  interval  of  N  periods  (during  which  N  components 
are  repaired)  is  described  by  the  binomial  distribution;  that 
is, 


p(x;N)  = 


i)p*  u-p)n-x 


(7) 


where  x=0,  1,  2,  .  ..,  N. 

However,  for  this  delivery  alternative  we  would  make 
a  delivery  only  if  at  least  one  demand  occurred  during  the  N 
periods.  Therefore,  we  must  condition  (7)  for  at  least  one 
demand  before  we  can  talk  about  the  expected  shipping  and 
shortage  costs.  Equation  (8)  provides  the  needed  form. 


p  ( x  ?  N ,  x  >1 )  = 


[n'i  ?*  d-p) 


[i  -  d-p) 


N-x 
IT 


;8) 


where  x  =  1 ,  2 ,  3 , 


N. 


Cost  Elements — As  with  the  deterministic  model  presented  in 
Chapter  I,  the  costs  which  will  be  considered  in  this  and  the 
other  delivery  alternatives  are  the  costs  associated  with 
delivering  the  item  to  the  NARF  and  the  costs  of  delay  asso¬ 
ciated  with  not  delivering  the  item  as  soon  as  it  is  demanded. 

Let  CT  represent  the  round-trip  costs  of  a  delivery. 
These  will  include  the  costs  of  the  truck  and  driver  from  the 
time  the  truck  starts  being  loaded  at  the  NSC  until  it  returns 
to  the  NSC.  CT  will  be  incurred  each  time  a  delivery  takes 
place.  If,  under  the  scheduled  delivery  scheme,  a  demand  does 
not  occur  during  an  interval  of  N  periods  then  we  assume 
initially  that  no  penalty  cost  tt  is  incurred  for  cancelling 
a  delivery.  Later,  we  will  incorporate  that  cost  into  the 
model . 

The  delay  costs  are  costs  incurred  at  the  NARF  as  a 
consequence  of  not  having  a  needed  part  on  hand  at  the  instant 
it  is  needed.  Two  elements  need  to  be  considered  here.  The 
first  is  a  cost,  denoted  by  S,  which  is  the  cost  associated 
with  putting  the  component  aside.  This  could  include  putting 
it  on  a  shelf  in  a  storage  area  and  documenting  how  far  repair 
had  progressed  and  what  had  been  ordered.  The  second  element, 
denoted  by  CD,  is  a  time  dependent  cost  which  we  will  refer 
to  as  the  delay  cost  per  period.  This  element  might  include 
labor  costs  due  to  work  stoppage,  inventory  holding  costs, 
and  costs  associated  with  a  repaired  component  not  being  avail¬ 
able  to  a  customer  such  as  a  fleet  squadron.  In  the  initial 
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formulation  of  the  expected  delay  costs  we  will  include  only 


the  CD  element. 

To  develop  the  expected  total  costs  per  period,  we 
must  first  develop  the  expected  total  costs  for  a  time  inter¬ 
val  of  N  periods,  given  that  at  least  one  demand  occurs  and 
hence  that  a  delivery  will  be  made.  We  have  already  identified 
CT  to  be  the  round-trip  delivery  cost.  We  must  now  determine 
the  expression  for  the  expected  total  delay  costs.  This  will 
be  the  product  of  CD  and  the  expected  total  delays. 

Expected  total  delays — The  delays  are  a  function  of  the  number 
of  different  configurations  that  demands  can  take  over  N 
periods.  For  example,  if  N  -  2  then  there  are  three  possible 
configurations.  The  one  where  a  demand  occurs  in  the  first 
period  and  none  occurs  in  the  second  results  in  a  delay  of 
one  period.  For  the  configuration  having  a  demand  in  the 
second  period  and  none  in  the  first,  a  delay  of  zero  periods 
results.  The  final  configuration  when  N  =  2  is  a  demand  in 
both  periods.  This  creates  a  delay  of  one  period  for  the 
first  period's  demand  and  zero  delay  for  the  second  period's 
demand,  resulting  in  a  total  delay  of  one  period. 

The  total  number  of  configurations  where  exactly  x 
demands  occur  is 


n 


x 


As  a  consequence,  the  total  number  of  configurations  which  can 
occur  over  N  periods  and  which  have  at  least  one  demand  is 
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n  =  Z  n 


=  2-1 


Theorem  II-l.  The  expected  total  delay  time  for  Alternative  1 


ETD(N)  = 


N(N-l)  p 
2[1' -  (l-p)N] 


Prop/.  To  determine  the  expected  total  delay  associated  with 
the  n  configurations  given  by  (9)  we  first  consider  only 
those  having  exactly  x  demands  where  x  >_  1.  The  probability 
of  each  such  configuration  is : 


P  (x,N)  = 


?x  (l-p)N~X 


1  ~  (l-p)‘ 


The  number  of  configurations  having  a  demand  in  period 


1  <  j  <  N  is 


m  -  fN-1' 

"  U-lj  ’ 


It  is  important  to  note  that  m  is  independent  of  j .  Those 
demands  occurring  in  period  j  will  have  to  wait  until  period 
N  for  delivery  and  hence  each  must  wait  N-j  periods.  The 
total  of  all  delays  for  those  configurations  having  x  demands 
is  therefore: 
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TD(x,N) 


N 

v 

j  =  l 


N— 1| 
[x-lj 


[N-j]  = 


f  N-l' 

lx-1. 


N  (N-l) 
2 - 


N 

l  (N-j) 

j-1 


N-ll 

[x-lj 


(13) 


Prom  (11)  and  (13)  it  easily  follows  that  the  total 
expected  delays  over  all  x  values  is  given  by  (14)  which  is 
identical  to  (10)  and  the  proof  is  complete. 

N 

ETD(N)  =  l  TD(x,N) P (x,N) 
x=l 

=  N (N-l) p  £ 

2[1  -  (l-p)N]  x^l 

N  (N-l)  p 
2[1  -  (l-p)N] 

Expected  Total  Costs  per  Period — In  developing  the  expected 
costs  per  period,  we  must  consider  a  cycle  which  can  include 
not  making  a  delivery  since  we  stated  earlier  that  Alternative  1 
has  a  scheduled  delivery  at  the  end  of  N  periods  if  there 
is  at  least  one  demand  during  this  time  interval. 

Theorem.  II-2.  The  expected  total  costs  per  period  for 
Alternative  1  ares 


N-l 
X— 1 


x-1 


(1-p) 


N-x 
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ECP(N) 


Prop/:  We  begin  with  time  zero  and  assume  that  at  least  one 

demand  occurs  in  the  first  N  periods.  The  total  expected 
costs  over  these  N  periods  is  then  made  up  of  the  delivery 
cost,  CT  ,  for  one  trip  and  the  total  expected  delay  costs 
which  is  the  product  of  CD  and  equation  (14) .  The  cost  per 
period  is  obtained  by  dividing  each  term  by  N.  The  result 
is  given  as  formula  (16) . 

CT  .  V**-1'? 

H  2 [i  -  a-p)*1: 

The  associated  probability  of  at  least  one  demand  during  the 

N 

first  N  periods  is  1  -  (1-p)  . 

Next  we  assume  no  demands  occur  during  the  first  N 
periods  and  at  least  one  occurs  during  the  next  N  periods. 
The  expected  cost  per  period  is  then: 

CT  CD(N-l)p 

+  4C1  -  (1-P)N]  * 

The  associated  probability  of  no  demands  in  the  first  N 
periods  and  at  least  one  in  the  second  N  periods  is 

(l-p)N  [1  -  (l-p)N]  . 
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The  general  formulas  for  no  demands  during  (k-l)N 


periods  and  at  least  one  during  the  last  N  periods  are  (17] 
and  (18) . 


CT  CD(N-l)p 

kN  2k [ 1  -  (l-p)Nj 


(1-P)(k-1)N  [1  -  (1-P)N1  • 


The  expected  costs  per  period  over  all  k  values  is 


therefore 


°°  ^T  Cp(N-l)p  (k-l)Nr  N-, 

ECP(N)  =  Z  =-1  +  - - - —  (l-p)(K  1,N[1  -  (l-p)Nj 

k=l  KN  2k[l  -  (1-p)  j 


C^Ll-d-p)  J  Cn(N-l)p 


Z  h  d-o)(k-1)N  • 
k=l  K 


The  summation  term  of  (19)  can  be  rewritten  as  follows: 


v  1  „  _4  (k-1)  N  1  r  1  ,, 

Z  r-  (1-p)  =  - K  “  j7  d-p) 

:=1  K  (1-p)  k=l  K 


1  :  ak 

"  a  k-1  k 


~  [-in(l-a)]  , 


where  a  =  (1-p)  .  The  series  sums  to  the  negative  of  Hn(l-a) 

since  a  <  1. 
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When  (20)  is  combined  with  (19)  the  result  is  (15) 
and  the  proof  is  complete. 

Determination  of  Optimal  N — Because  N  takes  on  only  integer 
values,  the  use  of  finite  differences  is  appropriate  for 
determining  optimal  N.  Optimal  N  is  that  value  of  N  which 
satisfies  the  following  relation: 

ECP(N-l)  >  ECP(N)  _<  ECP  (N+l)  . 

Equivalently,  optimal  N  is  the  largest  value  of  N  such 
that 

AECP(N)  =  ECP (N)  -  ECP(N-l)  <  0  . 

Unfortunately,  the  algebraic  expression  for  AECP(N)  which 
results  from  using  (15)  and  taking  the  difference  between 
ECP (N)  and  ECP(N-l)  is  just  as  complex  as  (15).  Therefore, 
a  numerical  evaluation  of  (15)  for  a  range  of  N  values 
appears  to  be  the  most  practical  way  of  searching  for  optimal 
N. 

Penalty  from  Cancelling  a  Scheduled  Delivery — Suppose  that 
the  cost  of  delivery  includes  a  penalty  cost  tt  for  cancel¬ 
ling  a  scheduled  delivery  (because  no  demands  occurred  during 
an  interval  of  N  periods) .  If  no  demands  occur  during  the 
first  (k-l)N  periods  then  k-1  cancellations  will  occur.  Since 
the  cost  for  each  is  tt  ,  equation  (17)  will  be  modified  to 
that  given  by  (21) . 
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(k-1)  7T  CD(N~I)  P 

kN  kN  2k[l  -  (l-p)N] 


The  expected  costs  per  period  over  all  k  values  is  then 


(C  -7T)Ei-U-p)N]  C  (N-1)P 

ECZP  (N)  =  „  +  -2 - * -  -  ~  2 - -  — !T  - 

(1-p) 


.  (22) 


The  impact  of  introducing  n  can  be  easily  seen  by  comparing 
(22)  with  (15). 

Delay  Cost  Independent  of  Time — Suppose  tjiat  there  is  a  fixed 
unit  cost  S  associated  with  putting  aside  a  component  until 
the  needed  part  arrives.  The  total  of  such  costs  during  an 
interval  of  N  periods  in  which  at  least  one  demand  occurs 
is  the  product  of  tr  and  the  expected  number  of  units  being 
delivered.  The  expected  number  of  units  being  delivered,  , 

is  given  by  (23) . 

„  _  ?  lx)  ?  (1'P} 

x  x=l  [1  -  (l-p)Mj 


[i  -  (i-P)N: 


[i  -  (i-p) 


Here  E(x)  is  the  expected  value  of  a  binomial  random  variable. 

The  total  expected  cost  over  N  is  then  increased  by 
the  product  SM^.  Formula  (17)  is  then  modified  as  follows 
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Sp  +  CD(N-l)p 


(24) 


kN  2k[ 1  -  (l-p)N] 


and  (15)  is  therefore  modified  as  shown  by  (25) 


ECP(N)  = 


C  El-(l-p) N]  C_ (N— 1) p 


N 


s-f 


-£n[l-(l-p)N] 


(1-P) 


N 


(25) 


If  the  time-dependent  delay  cost  CD  is  negligible 
relative  to  S  then  optimal  N  will  be  infinite  since  the 
S  term  in  the  left  bracket  of  (24)  is  not  an  increasing 
function  of  N. 

Time  Standards — Time  standards  have  been  established  as  upper 
bounds  on  the  average  time  allowed  a  supply  center  to  respond 
to  a  customer's  demand.  If  T  denotes  this  standard  in 
periods  then  the  expected  delay  in  periods  per  unit  delivery 
is  required  to  not  exceed  T. 

Theorern  II-3.  The  expected  delay  per  unit  under  alternative  1 
is 

UD(N)  =  .  (26) 


Proo/.  Equation  (13)  gives  the  total  of  all  delays  for  those 
configurations  having  x  demands  during  N  periods. 

Equation  (11)  gives  the  probability  of  each  such  configuration. 
Dividing  (13)  by  x  gives  the  total  of  all  average  delays  per 
unit  for  those  configurations.  Thus  the  expression  for 
UD(N)  is 
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UD(N) 


P (x,N) 


N 

Z 

X=1 


TD(x,N) 

x 


N(N-l) p  J  1 

2 [1  -  (l-p)N]  x=l  x 


rN-ll 

[x-lj 


P*’1 


(1-p) 


N-x 


(N-l) 

2 [ 1  -  (l-p)N] 


N 

I 

x—  1 


Uj 


x , .  , 

p  (1-p) 


N-x 


(N-l) 

2 [ 1  -  (l-p)Nj 


[1 


d-p) 


N, 

J 


Cancellation  of  the  probability  terms  in  the  denominator  and 
numerator  results  in  the  right  side  of  (26)  and  the  proof  is 
complete. 

Setting  up  the  time  constraint  inequality;  namely, 


N-l 

UD  (N)  =  <_  T  ; 


and  solving  for  N  gives 


N  <  2T  +  1  (27) 

Inequality  (27)  provides  an  easy  way  of  determining  if  an 
optimal  N  is  also  feasible  with  respect  to  a  given  time 
standard. 

Truck  Capacity  Constraint — When  CT  was  introduced  it  was 
assumed  to  represent  the  round-trip  delivery  costs  incurred 
by  one  truck  (and  driver) .  As  a  consequence,  it  is  appropriate 
to  include  a  constraint  to  represent  the  capacity  of  that 


4 
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truck.  Let  M  represent  the  maximum  number  of  units  of  a 
given  repair  part  that  can  be  loaded  on  a  given  truck.  Under 
alternative  1  it  is  possible  to  have  a  demand  during  each  of 
a  sequence  of  N  periods.  As  a  consequence,  N  units  would 
need  to  be  delivered.  Therefore,  N  <_  M  is  required  to  insure 
that  all  possible  quantities  demanded  during  N  periods  can 
be  delivered  by  the  truck  at  the  end  of  these  periods. 
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III.  ALTERNATIVE  2 


Introduction — Alternative  2  assumes  that  a  delivery  is  not 
made  until  K  units  of  an  item  have  been  demanded.  Delivery 
is  assumed  to  take  place  as  soon  as  the  last  demand  occurs. 
The  basic  assumptions  spelled  out  in  Chapter  II  with  respect 
to  cost  elements  and  the  probability  of  demand  for  a  repair 
part  also  apply  to  this  alternative. 

Probabilities — For  Alternative  2  we  will  be  interested  in  the 
probability  that  n  periods  will  pass  before  we  accumulate 
K  demands  with  the  last  demand  occurring  in  period  n.  This 
probability  is  described  by  the  negative  binomial  distribu¬ 
tion  [  2]  ?  namely. 


p(n;  K) 


n-1 ' 

(k-iJ 


K 


P 


(1-p) 


n-K 


(28) 


where  n=K,  K+l,  K+2 ,  ...  . 

Expected  Total  Delays — As  was  observed  in  Chapter  II,  the 
delays  are  a  function  of  the  number  of  different  conf igurations 
that  demands  can  take  over  n  periods .  When  K  =  2  and  n  =  4 , 
we  get  3  demand  configurations.  In  all  three,  the  second 
demand  occurs  in  the  fourth  period.  Thus,  the  first  demand 
can  occur  in  either  period  1,  2 ,  or  3 .  The  delay  is  zero  for 
the  second  demand  in  all  three  cases.  The  delay  for  the  first 
demand  is  3,  2,  and  1  periods,  respectively. 

For  the  general  case,  the  number  of  configurations, 
m,  having  a  demand  occurring  in  periods  1  <  j  <  n  -  1  is  given 
by  (29). 
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=  fn-a 

(k-2 


Equation  (29)  differs  from  (12)  because  the  K  demand  always 
occurs  in  period  n  for  Alternative  2.  Again  we  see  that  m 
is  independent  of  j . 

Theorem  111*1.  The  expected  total  delay  time  for  Alternative  2 


ETD(K)  =  , 


Proof .  The  total  delay  associated  with  the  conf igurations 
having  a  demand  in  period  j  is  m[n-j].  The  total  of  all 
delays  for  K  demands  during  n  periods  is  therefore 


TD(n,  K)  =  I  (n-j)fc? 

3-1  lK  2 


We  see  that  (31)  can  be  rewritten  as 


TD(n,  K)  =  [S“?]  E  (n-j) 
l  J  -i=l 


n(n-l)  f  n—  2 1 

2  (K-2j 


n(K-l)  fn-1 
2  K-l 


The  probability  of  each  configuration  is 


P(n,K)  =  pK(l-p)n"K 


20 


The  total  expected  delays  over  all  n  values  are 
obtained  by  summing  the  products  of  (32)  and  (33)  over  all 
n  >  K. 


ETD(K)  =  Z  TD (n,K) P (n,K) 
n=K 


n=K 


n(K-l)  f n-l] 
2  [k-1 


X  . .  , n-K 

p  (1-p) 


OO  f 

K-1  „  fn-1 

ni/l''-1 


K..  ,  n-K 

P  (1-p) 


K-1  “  r  ^  f K+x-ll  K  ..  , ; 

^0C  *  K-  l  k-1  JP  (i-P' 


K-1 

2 


K-1 

2 


K  + 

r 

Ik  + 


:  f K+x-i 

"aX1  k-i 
x=0 

K  ( 1-p) 


K ,,  , 

p  (1-p) 

(K-1) 

2p 


x 


(34) 


Since  (34)  is  identical  to  (30),  the  proof  is  complete. 
Expected  Total  Costs  per  Period — Under  Alternative  2  a 
delivery  will  be  made  as  soon  as  K  demands  have  accumulated. 
This  will  span  n  ^  K  periods. 

Theorem.  III-2.  The  expected  total  costs  per  period  for 
Alternative  2  are : 


ECP(K) 


00 


l 

n=K 


1  [n-l] 
n |K-lJ 


K . 
p  (1-p) 


n-K 


CD(K-1) 

2 


(35) 
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Proof .  For  a  specified  value  of  K,  the  delivery  costs  per 
period  are 


ETCP (K) 


n=K 


n-l) 

[K-1JP 


K,,  .  n-K 

(1-P) 


(36) 


Equation  (36)  results  from  summing  the  product  of  C^/n  and 
(28)  over  all  n  ^  K. 

The  expected  delays  costs  per  period  are  obtained  by 
first  dividing  the  total  delays  for  a  given  n  (namely 
TD(n,K))  by  n  and  then  summing  the  product  of  CD  ,  that 
result,  and  (33)  over  all  n  >  K. 


EDCP(K)  =  C  Z  p  (n,K) 

n=K  n 


=  C 


D 


(K-l) 

1 


oo 

z 

n=K 


"n-ll 

IK-lJ 


K  , .  , 

P  (1-p) 


n-K 


=  C 


D 


(K-l) 

2 


(37) 


Adding  (36)  and  (37)  together  gives  (35) . 

Determination  of  Optimal  K — The  use  of  finite  differences  to 
determine  optimal  K  does  not  result  in  any  relationship 
which  is  less  complex  than  (35) .  In  particular,  the  infinite 
sum  in  the  term  remains.  Fortunately,  however,  bounds 

can  be  derived  for  (35)  which  can  serve  to  reduce  the  number 
of  K  values  for  which  (35)  must  be  evaluated  in  searching 
for  optimal  K. 
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Tkzaxtm  III-3.  The  expected  delivery  cost  per  period  for 
Alternative  2  is  bounded  as  follows: 


CTP  C  p 

-±-  <  ETCP  { K )  <  .  (38) 


P-ro  o  £ .  Equation  (36)  is,  in  reality,  the  product  of  CT  and 
the  expected  value  of  1/n;  that  is, 


Jensen's  Rule  [2]  therefore  allows  us  to  obtain  the  lower 
bound.  That  rules  states  that 


ETnT 


Because  CT  is  a  constant  we  can  write 

("1>  (  i ) 

e TKT  ±  cte(HJ  ‘  ETCP<K)  ■ 

The  value  of  E(n),  when  n  has  the  negative 
binomial  distribution  given  by  (28),  is,  from  [2], 


E(n,  - § . 


(39) 


Therefore, 


cT  cTp 

ETnT  =  ~K~ 


(40) 
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Equation  (40)  is  the  lower  bound  of  (38) • 
Next, 


Z 

n=K 


rn-l] 

K,.  .  n-K  ”  n-1 

fn-21 

IK-IJ 

P  (1-P)  =  2 

n=K 

Ik-i 

K,,  x 
P  (1-p) 


n-K 


<  Z 

n=K 


n-2 

K-l 


K/n  v  n-k 
p  (1-p) 


since 


n-1 

n 


< 


n 

n 


1 


Then, 


00 

z 

n-K 


fn-2] 

K-l 

\  y 


K,,  x 
P  (1-p) 


n-K 


PK-iu-pi 


n-K  _  p 
~  K-l 


Therefore, 


E 


T 

n. 


and  hence  the  upper  bound  of  (38)  is  (41) 


(41) 


When  (38)  is  combined  with  (37) ,  the  bounds  for  the 
expected  total  costs  per  period  are  easily  determined  to  be: 


CTP  Cn[K-l]  C  p  C  [K-l] 

—  -2-2 -  <  ECP(K)  <  jSj.  +  -=~2 - 


(42) 
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These  bounds  can  be  analyzed  using  finite  differences  to  deter¬ 
mine  their  respective  optimal  K  values.  Theorem  III- 4  and 
its  proof  give  the  results. 

‘Theorem  III-4.  The  value  of  K  which  minimizes  the  lower 
bound  of  (42)  is  one  less  in  value  than  the  value  of  K  which 
minimizes  the  upper  bound. 

Proof .  Using  the  approach  of  finite  differences  we  know  that 
optimal  K  for  the  lower  bound  is  the  largest  K  value  for 
which 

C  (K)  -  C(K-l)  <  0  . 

When  we  evaluate  this  difference  we  get 

CTP  CD 

C(K)  -  C(K-l)  =  -  xjjfZTT  +  ~T  <  0  ' 
which  can  be  rewritten  as 

2C  p 

K(K-l)  <  .  (43) 

Similarly,  for  the  upper  bound 

Gw  -  c(K-i)  =  -  ~(K-iy~(K-2)  +  -r  <  0  ' 

which  can  be  rewritten  as 

2C  P 

(K— 1)  (K-2)  <  -jJ-  .  (44) 
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We  see  immediately  that  (43)  and  (44)  differ  by  the 
fact  that  K  has  been  replaced  by  K-l.  As  a  consequence, 
optimal  K  for  (43)  will  be  one  unit  less  than  optimal  K 


for  (44)  . 

Conc_eature^ .  The  value  of  K  which  minimizes  ECP(K)  is  either 
K*  or  K**  where  K*  minimizes  the  lower  bound  of  (42)  and 
K**  minimizes  the  upper  bound. 

Although  computational  experience  shows  this  conjec¬ 
ture  to  be  true  and  intuition  strongly  favors  it,  no  formal 
proof  has  been  yet  been  discovered. 

Modification  of  the  Expected  Costs  per  Period — In  Chapter  II, 
two  modifications  of  the  basic  expected  costs  formula  were 
presented.  The  first  was  the  inclusion  of  a  penalty  cost  for 
cancelling  a  scheduled  delivery.  Such  a  cost  would  not  be 
appropriate  to  Alternative  2  since  it  is  an  unscheduled  de¬ 
livery  strategy.  However,  it  could  be  argued  that  C T  for 
this  alternative  might  be  appropriately  higher  than  that  for 
a  scheduled  delivery  since  the  need  for  a  truck  will  not  be 
known  until  the  Kth  demand  has  occurred. 

The  second  modification  was  the  inclusion  of  a  fixed 
delay  cost  S  per  unit  demanded.  Under  Alternative  2,  the 
total  of  such  costs  is  merely  SK  since  K  units  are  always 
delivered.  The  expected  costs  per  period  is  then  modified  to 
include  the  term 


SK 


Z 

n=K 


lfn-r 

n[K-lJ 


K 


P 


(1-p) 


n-K 
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Thus  (35)  is  now 


G°  ,  . 

ECP(K)  =  [C_  +  SK]  T.  i(n_1l 


‘J  n=K  n|K-1 


Kf.  .  n-K  Cd(K_1) 
P  (1-P)  +  - * - 


(45) 


The  bounds  for  ECP(K)  are  modified  also.  The  lower 


is  now 


C  p  C  (K-l) 

IT  +  SP  +  - J— 


and  the  upper  is 


CTP  ,  SpK  J  CD(K_1) 

K=T  K— 1  2 

If  the  value  of  CQ  is  negligible  with  respect  to  S  then 
the  value  of  K  which  minimizes  the  lower  bound  is  infinity. 
The  value  of  K  which  minimizes  the  upper  bound  is  also 
infinity.  Therefore,  optimal  K  minimizing  ECP(K)  must  also 
be  infinite. 

Time  Standards — Again  denote  T  as  the  upper  bound  on  the 
average  delay  time  per  unit. 

Theorern  III-5.  The  expected  delay  per  unit  under  Alternative 
2  is 

UD  (K)  =  .  (46) 

Prop/.  Because  K  units  will  always  be  delivered  we  can  get 
the  expected  delay  per  unit  by  merely  dividing  the  total  ex¬ 
pected  delay,  given  by  (30),  by  K.  The  result  is  (46). 
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When  the  time  constraint  is  introduced,  we  get 


UD(K)  =  <  T  . 

Solving  for  K  in  this  inequality  gives 

K  <_  2pT  +  1  .  (4  7) 

Truck  Capacity  Constraint- — If  M  is  the  maximum  capacity  of 
the  truck  then  it  follows  that  K  <  M  is  the  capacity  constraint. 
Expected  Number  of  Periods  Between  Deliveries — When  comparing 
Alternative  2  with  the  other  alternatives  it  will  be  useful  to  know 
the  expected  number  of  periods  between  deliveries.  The  formula 
was  given  earlier  as  (39);  namely, 

E(n)  =  |  . 
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IV.  ALTERNATIVE  3 


Introduction — Under  Alternative  3,  we  start  counting  time  from 
when  the  first  demand  occurs  after  the  delivery  truck  has 
returned  and  is  again  available  for  further  deliveries.  We 
then  wait  N-l  periods  after  the  first  demand  before  we  deliver 
again.  The  basic  assumptions  spelled  out  in  Chapter  II  rela¬ 
tive  to  cost  elements  and  the  probability  of  demand  for  a 
repair  part  also  apply  to  this  alternative. 

Probabilities — The  probability  of  x  demands  in  N  periods, 
given  that  the  first  one  always  occurs  in  the  first  period,  is 


P  (x;N) 


N-l1 

.X-l, 


(48) 


where  x  =  1,  2,  3, 

Expected  Total  Delays — We  start  counting  delay  time  from  the 
period  when  the  first  demand  occurs.  If  we  deliver  N-l  periods 
later  then  the  total  delay  for  the  first  demand  is  N-l  periods. 
The  delays  associated  with  subsequent  demands  occurring  in 
periods  2,  3,  up  to  N,  are  a  function  of  the  different  con¬ 
figurations  of  demand  which  can  occur.  For  example,  if  N=2 
then  two  cases  occur;  the  first  has  a  demand  in  period  1  and 
none  in  period  2,  the  second  has  a  demand  in  both  periods  1 
and  2.  In  each  case  the  first  demand  is  delayed  one  period. 

In  the  first  case,  no  subsequent  delay  occurs  since  no  sub¬ 
sequent  demand  occurs.  In  the  second  case,  the  second  demand 
incurs  a  delay  of  essentially  zero  since  the  truck  leaves 
shortly  after  the  demand  occurs. 
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The  total  number  of  configurations  having  exactly  x 
demands  during  N  is 


n 


x 


N— 1^ 
x-1. 


(49) 


The  total  number  of  configurations  for  a  given  N  is  then 


n 


N 

Z 

x=l 


fN-ll 

x-1 


2 


N-l 


(50) 


Theorem.  IV- 1.  The  expected  total  delay  time  for  Alternative  3 
is : 


ETD(N)  =  [N-l] 


(N-2 ) 


p  +  1 


(51) 


Prop/.  The  number  of  configurations  having  a  demand  in  period 
2  <  j  <  N  and  a  total  of  x  demands  over  the  N  periods  is 


m  = 


N-2] 
x-2  j 


(52) 


where  x  =  2,  3,  .  .  . ,  N.  Obviously,  m  is  independent  of  j. 
The  delay  for  a  demand  occurring  in  period  j  is  N- j . 

The  delay  for  demands  occurring  in  the  first  period 
is  N-l  and  equation  (49)  gives  the  total  number  of  configura¬ 
tions  having  x  demands.  The  total  of  all  delays  for  these 
x  demands  is 
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TD(x,N) 


1 1  X  [ 

[x-lj 

(N-l) 

+  S 

j=2  1 

fN-1] 

[  x- 1 J 

(N-l) 

fN-2  j 
(x-2  j 

f  N-l] 
[x-lj 

(N-l) 

+ 

& 

1  1 

M  h-* 

N  (N 

2 

_ 1 

f  N-l"} 

[x-lj  * 

in  —  x  ; 

x-2  J 


(N-j) 


(N-l) (N-2) 


(x-1) (N-2) 
2 


(53) 


The  probability  of  one  configuration  having  x  demands 
in  N  periods,  given  that  one  occurs  in  the  first  period  is 


_  ,  ...  x-1  ,,  .N-x 

P (x,N)  =  p  (1-p) 


(54) 


The  expected  total  delay  is  then  obtained  by  summing 
the  products  of  (53)  and  (54)  over  all  1  <  x  <  N. 

N 

ETD(N)  =  Z  TD (x , N)  P(X,N) 
x=l 


_  r  Tn  .  (N-2)xl 

-  x-1  U 


N-l!  x-1,,  . N-x 

x-lj  P  (1-P} 


N  N 
—  z 

2  i 

X=1 


N-l 

x-1 


x-1,.,  ,  N-x 

p  d-p) 


N-2  N 
+  Z  x 


x=l 


px-1(i-p)'1-x 
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N  .  N-2 
=  I  +  ~ 


N-l 

Z 

u=0 


(u+1) 


N-l' 
u  j 


pu(l_p)N-i-u 


N  .  N-2 

=  2  +  ~T 


(N-l) 


p  +  1 


=  (N 


-1)  1  + 


N-2 


The  proof  is  complete. 

Expected  Total  Costs  per  Period — In  analyzing  the  expected 
costs  per  period  we  consider  a  cycle  whose  length  is  N  plus 
the  time  span  in  periods  between  the  time  the  truck  becomes 
available  for  a  delivery  and  the  first  demand  occurs. 

Theorem  IV-2.  The  expected  total  costs  per  period  for 
Alternative  3  are 


ECP(N)  =  ETC (N) 


+N  P(1“P} 


k-1 


where : 


ETC (N)  =  CT  +  Cd(N-1) 


(N-2)  n  +  1 
— 2 -  P  +  1 


(55) 


(56) 


Prop/.  Equation  (56)  describes  the  expected  total  costs 
over  the  N  periods.  Now,  if  the  first  demand  occurs  in  the 
first  period  after  the  truck  becomes  available  then  the  total 
cost  per  period  is  ETC(N)  divided  by  N  since  we  deliver  in 
the  Nth  period.  The  probability  of  a  demand  in  the  first 
period  is  p.  If,  instead,  the  first  demand  occurs  in  period  2, 
then  the  average  expected  cost  per  period  is  ETC(N)  divided 
by  N+l  and  the  probability  of  the  first  demand  occurring  in 
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* 


t±' 


1 


period  2  is  (l-p)p.  The  general  forms  for  the  average  cost 
per  period  and  the  associated  probability  when  the  first  demand 
occurs  in  period  k  are: 


ETC (N) 
(k-1) +N 


and 


M  .  k-1 
p(l-p) 


The  expected  total  costs  per  period  are  the  sum  of  the  products 
of  these  two  general  forms  over  all  possible  k  values.  The 
result  is  equation  (55) .  The  summation  term  of  equation  (55) 
cannot  be  written  in  a  closed  form,  however,  bounds  can  be 
stated. 

Theorem  IV-3.  The  expected  total  costs  per  period  under 
Alternative  3  are  bounded  as  follows: 


ETC (N) 


irrfhjp 


<  ECP(N)  <  ETC (N) 


min 


lN 


(-*np) 


I 

*J 


(56) 


Prop/.  We  can  write  ECP(N)  as  the  product  of  ETC (N)  and  the 
expected  value  of  the  reciprocal  of  (N+k-1)  where  k  has  a 
geometric  distribution.  Now  from  Jensen’s  Rule  [2]  we  know 
that 


1 

eTn+FTT 


Next, 


E (N+k-1)  =  N— 1  +  E (k) 

1  +  (N-l)p 
P 


(57) 


(58) 


33 


A 


The  product  of  the  lower  bound  from  Jensen's  Rule  and  ETC(N) 
is  the  lower  bound  given  by  (56);  namely, 

ETcm  lVfg-Up  • 

Consider  next  the  actual  expansion  of  the  bracket 
term  of  (55) . 

E  N+k-1  =  N  P  +  N+I  P  (i_P)  +  N+2  P  (1_P)  +  * 

When  N=1  the  right  side  of  (60)  reduces  to 

p  +  P<1-P».  *  p'lrp)?  +  pg-p).3  ♦  — 


=  I  [-£n (1-a) ] 

cl 

=  itp  C“£np]  . 

Here  we  introduced  a  in  place  of  (1-p)  to  simplify  the 
middle  two  steps.  When  N=2  we  get 

j  P  +  j  P(l-P)  +  \  P(l-P)2  +  — 


=  [-  £n(l-a)  -  a]  . 

& 


(59) 


(60) 


(61) 
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For  general  N: 


E  f - - — ^  =  £— 

(N+k-lj  aN 


N-l  k 
tn(l-a)  -  Z  3_ 

k=l  *  J 


(62) 


Equation  (62)  can  be  used  in  place  of  the  infinite  sum  for 
computations  of  ECP(N). 

An  upper  bound  is  suggested  by  (62) ;  namely, 


P 

(1-P) 


N 


L-inp]  . 


N 

However,  (61)  provides  a  tighter  bound  since  1  -  p  >  (1-p) 
for  N  >_  2. 

Inspection  of  (60)  suggests  another  upper  bound.  If 
we  keep  all  denominators  at  a  value  of  N  we  get 

|  fl  +  (1-p)  +  (1-P)2  +  (1-P)3  +  - ] 


(63) 


The  advantage  of  (63)  over  (61)  is  that  it  is  a  de¬ 
creasing  function  of  N  whereas  -2,np  remains  constant.  When 
p  <  0.368,  (63)  gives  the  best  result.  However,  when  p 

increases  to  0.7  or  0.9,  we  need  N  values  of  at  least  3  and 
10,  respectively  before  (63)  is  preferred.  Thus,  we  can  state, 
in  general,  that  our  bound  should  be 

min  ,  -£np}  (64) 
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And,  when  we  form  the  product  of  (64)  with  ETC(N)  we  get 


ETC (N) 


I^p  min 


(-inp)  } 


which  is  the  upper  bound  of  (56) . 

While  (62)  is  not  difficult  to  use  when  N  is  small, 
the  summation  becomes  tedious  for  large  N  and  hence  optimal 
solutions  to  the  bounds  should  be  valuable  in  narrowing  the 
values  of  N  to  be  investigated.  As  in  Chapter  III,  we  have 
not  been  able  to  prove  it  but  computational  experience  suggests 
that 

•  The  value  of  N  which  minimizes  ECP(N)  is  bounded 
as  follows: 


N*  <_  N  <_  N**  , 

where  N*  minimizes  the  lower  bound  and  N**  minimizes  the 
upper  bound  given  in  (56) . 

The  relationships  for  determining  the  optimal  values 
of  N  for  the  lower  and  upper  bounds  are  derived  below. 

Theopern  IV-4.  The  optimal  value  of  N  which  minimizes  the 
lower  bound  is  the  largest  value  of  N  which  satisfies 

(N-l)  (N-2)  pZ  +  2[  (N-2)  p  +  1]  <  .  (65) 


The  optimal  value  of  N  which  minimizes  the  upper  bound  is 
the  largest  value  of  N  which  satisfies 
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N (N-l) p  +  2 ( 1-p) 


(66) 


Proof.  For  the  lower  bound  we  want  the  largest  N  such  that 


ETC(N)  T+Tn-TTp  '  BTC(H-1) 


[l+(N-2)p] 


The  difference  in  the  CT  term  is 


(N-2) p] 


~CTP 


The  difference  in  the  CD  term  is 


p(N-l)[Sj^  p+l]  p(N-2)[Sjl  p+l 


D  1+  (N-l)  p 


1+ (N-2) p 


(N-l) (N-2 ) p  +  2 (N-2) p  +  2 
CDP  L1+(N-1) pjLl+N-2) pj 


The  sum  of  (68)  and  (69)  constitutes  the  left  side  of 
(67) .  Upon  cancellation  of  common  terms  and  movement  of  the 
cost  parameters  to  the  right  side,  we  get  the  result  shown  as 


For  the  upper  bound  we  assume  that 


min  (i  ,  -  fcnp}  =  i 
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for  the  N  values  where  the  conjecture  would  be  useful.  We 
then  want  the  largest  value  of  N  for  which 


ETC(tl)  [nt^pt]  -  etc(n-i)  [THrrfn=pr]<  0  •  (70) 

The  difference  in  the  CT  term  is 


ctp  ri  _  11.  'ctp 

(l-p)  Ln  N-l  J  N(N-l)(l-p) 
The  difference  in  the  CD  term  is 


[(N"2)P  +  2^  -  |£rf  C  (N-3)  p  +  2]] 


2N (N-l) (l-p) 


[pN (N-l)  +  2  (l-p)] 


The  sum  of  (71)  and  (72)  constitutes  the  left  side  of  (70)  . 
Cancellations  and  rearrangements  yield  (66) . 

Modification  of  the  Expected  Costs  per  Period — The  inclusion 
of  the  penalty  cost  for  cancelling  a  scheduled  delivery  is 
not  appropriate  for  this  alternative  since  it  is  not  a  sched¬ 
uled  delivery  strategy.  As  in  Chapter  III,  it  might  be  argued 
that  the  CT  value  would  be  higher  for  this  alternative  than 
for  Alternative  1.  It  might  also  be  argued  that  the  CT 
value  is  lower  than  for  Alternative  2  since  the  time  at  which 
a  truck  will  be  needed  is  known  as  soon  as  the  first  demand 
occurs  after  the  truck  completes  the  previous  delivery.  Thus 
a  reservation  can  be  made  ahead  of  time. 
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The  second  modification  was  the  inclusion  of  the 


fixed  cost  S  per  unit  demanded.  As  with  Alternative  1, 
the  total  expected  costs  ETC (N)  over  N  periods  is  the 
product  of  S  and  the  expected  number  delivered.  And  for 
Alternative  3,  the  expected  number  delivered,  ,  is: 


m3  = 


N 


x=l 


(N- ll _x-l 
x-1 


(1-p) 


N-x 


-  V  Lu.li"-1' 

u=0 


u  . N-l-u 

p  (1-p) 


=  1  +  E(U) 


=  1  +  (N-l)p 


(73) 


The  modified  form  of  ETC(N)  is  given  by  equation  (74) . 


ETC (N)  =  CT  +  S  +  (N— 1 ) 


CD(N-2)p 


Sp  +  Cr 


(74) 


The  lower  bound  for  ECP(N)  when  the  S  term  is 
included  is: 


pCT  +  pCD(N-l) 

♦  i] 

i  +  (N-: 

UP 

(75) 


and,  if  the  value  of  CQ  is  negligible  with  respect  to  S, 
the  optimal  value  of  N  for  this  bound  is  infinite.  A  similar 
argument  results  in  an  infinite  optimal  N  for  the  upper 
bound.  As  a  consequence,  the  N  minimizing  ECP(N)  must  also 


be  infinite.  The  reader  will  recall  that  this  was  also  the 
result  for  Alternative  1. 

Time  Standards — With  T  being  the  time  constraint  in  periods 
we  again  must  have  the  expected  delay  per  unit  such  that 

UD  (N)  <_  T  , 

where  UD(N)  is  the  expected  delay  per  unit. 

Theorem  IV-5.  The  expected  delay  per  unit  under  Alternative 
3  is 


UD (N)  = 


(n— 2 )  +  LkiLkiBl N '2 
p 


(76) 


Proof .  Following  the  arguments  in  the  proof  of  Theorem  II-3, 


UD (N)  =  E  TD.(.?'*ii.  P(x,N) 
x=l 


_  r  Tn  .  (N— 2 )  1  f N-l)  _x-l , ,  _ i  N- 

~  xJl  |_23?  +  -T-Jlx-lJP  (1'P) 


N 


1)  .  N-2 


(77) 


1 [ N— 1 I x— 1 . .  .N-x 
x|x-l  P  (1'P) 


(1-p) 


N-x 


_  i-(i-p) 

p 


Substitution  of  (78)  into  (77)  results  immediately  in  (76) . 

In  contrast  to  Alternative  1,  substitution  of  (76) 
for  UD(N)  in  the  time  constraint  does  not  provide  a  simple 
expression  for  N. 

Truck  Capacity  Constraint — If  M  is  the  maximum  capacity  of 
the  delivery  truck  then  the  truck  capacity  constraint  is 
N  <_  M  since  it  is  possible  that  N  units  will  be  demanded 
over  the  N  periods . 


V.  ON-SITE  SYSTEM 

Introduction — An  on-site  inventory  system  is  an  alternative 
to  delivery.  In  comparing  it  to  direct  delivery,  it  is  appro¬ 
priate  to  develop  the  comparable  total  expected  costs  per 
periods.  Therefore,  an  expected  costs  model  for  the  on-site 
system  must  be  developed  for  comparison  with  each  of  the  three 
direct  delivery  alternatives. 

A  basic  assumption  of  these  models  will  be  that  suf¬ 
ficient  inventory  exists  on-site  to  meet  all  demands. 

Cost  Elements — A  fixed  cost  per  period  is  assumed  to  represent 
labor  and  related  overhead  costs.  It  will  be  denoted  by  . 

A  delay  cost  is  also  appropriate  since  it  is  assumed 
that  the  on-site  system  will  not  be  distributed  to  each  re¬ 
pair  production  line  but  will  instead  be  centrally  located  at 
the  NARF.  This  delay  cost  is  assumed  to  be  a  fixed  cost  per 
unit  demanded  and  will  be  denoted  by  s. 

Alternative  1  Model — In  comparing  the  expected  costs  of  the 
on-site  system  with  the  scheduled  delivery  alternative,  we 
need  to  evaluate  the  expected  delay  costs  over  the  same  time 
frames  as  were  used  in  Alternative  1.  If  we  have  k-1  inter¬ 
vals  of  N  periods  before  a  demand  occurs  then  the  probabil¬ 
ity  of  such  an  occurrence  is  given  by  formula  (18) .  The  ex¬ 
pected  number  of  units  demanded  during  the  last  interval  of 
N  periods  is  given  by  equation  (23)  and  the  expected  total 
delay  costs  are 
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(79) 


SM,  =  - ^ - rj- 

X  r-«  .  -  \N-» 

Ll  -  (1-p)  J 


The  expected  delay  costs  per  period  over  the  kN  periods  is 
then 


kN  .  ^  v  N-j 

kLl  -  (1-p)  ] 


The  expected  value  of  (80)  over  all  possible  k  values  is 
determined  from 


Cl  -  (1-P)N]  k-l  K 


£  p(l-P)  (k~1)Nn  -  (i-p)N] 


„  L  1 ,,  ,  (k-l)N 

=  sp  k-i  £(1*p) 


-Mi  -  (i-P)  l 


(1-p)' 


Finally,  the  total  expected  costs  per  period  for  the 
on-site  system,  OECP(N),  is  the  sum  of  CL  and  (81);  namely. 


OECP(N)  =  CL  +  sp  -  Jill 

(1-p) 


Alternative  2  Model — The  CL  term  is  again  appropriate.  The 
total  delay  cost  is  now  the  product  sK  for  one  delivery 
cycle.  The  delay  costs  per  period,  where  n  is  the  number 
of  periods  required  to  accumulate  K  demands ,  is 
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The  probability  of  n  periods  being  required  is 


given  by  (28)  . 

Combining  (28)  and  (83)  and  summing  over  all  possible 
n  values  gives 


SK 


Z 

n=K 


1  fn-11 
n^K-lj 


PK(1-P) 


n-k 


(84) 


And  the  total  expected  costs  per  period  is  the  sum  of  CT 

Li 

and  (84);  that  is, 


OECP(K) 


CL  + 


SK 


Z  i 
n=K 


n-1 

K-l 


PK(1-P)n~k 


(85) 


Alternative  3  Model — The  expected  number  of  units  demanded 
before  delivery  under  Alternative  3  is  given  by  (73).  Therefore, 
the  total  delay  costs  are 


s[l  +  (N-1) p] 


The  number  of  periods  between  deliveries  is  (k-l)  +  N  so  that 
the  delay  cost  per  period  is 


s[l  +  (N-1) p] 
- (k-l)  +  N  ' 


(86) 


and  the  probability  of  the  first  demand  occurring  in  period 
k-l 

k  is  p(l-p)  .  The  expected  value  of  (86)  over  all  possible 


k  values  is  therefore 


VI .  TRADE-OFF  MODELS 

Introduction — The  first  trade-off  models  to  be  presented 
below  compare  each  direct  delivery  alternative  with  its  asso¬ 
ciated  on-site  model.  These  models  include  the  special  con¬ 
ditions  where  optimal  N  and  K  are  unity  since 
they  appear  to  correspond  to  situations  under  which  an  on¬ 
site  system  seems  most  appropriate  (i.e.,  having  an  on-site 
system  seems  better  intuitively  than  making  a  direct  delivery 
every  time  a  demand  occurs) .  Finally,  the  models  incorporate 
the  time  standard  and  truck  capacity  constraints. 

The  individual  models  are  then  combined  into  a  com¬ 
posite  model  which  is  designed  to  resolve  the  question  of 
which  direct  delivery  alternative  to  use  if  direct  delivery 
is  the  optimal  strategy. 

Alternative  1 — When  comparing  the  direct  delivery  costs  given 
by  (15)  with  the  corresponding  on-site  costs  given  by  (82) , 
it  follows  that  the  on-site  system  is  preferred  when 
OECP(N)  <  ECP(N).  When  the  cost  elements  are  introduced  into 
that  inequality,  the  result  is  that 


(88) 

Now  (15)  is  minimized  by  a  certain  N  value.  However, 
the  right  side  of  (88)  is  minimized  for  a  slightly  larger 
value  of  N  because  of  the  (-2s)  term  creating  a  savings  in 
the  delay  costs.  It  is  this  latter  value  that  is  appropriate 


<  CT[l-(l-p)N]  |  rC„(N-l)  -2s]p 


(1-P) 


when  considering  the  trade-off  between  on-site  and  direct 
delivery. 

To  have  N  =  1  minimize  the  right  side  of  (88)  we 
need  to  have  the  latter's  value  be  less  for  N  -  1  than  for 
N  =  2.  When  N  =  1,  (88)  reduces  to 


(39) 


When  N  =  2,  the  right  side  of  (88)  is  reduced  as  fellows: 


T 


:i-(i~P)2: 


-  rvl-  ( l-p)  2 : 

a-?)2 


c  p  c  p 

~2~  (2-p)  +  -J-  - 

sp 

- Inp ( 2-p) 

- 

- 

l  (l-p) 

CTp  CTP  CDP 

T  -  SP  -  2  +  2 

-£np(2-p) 

-  d-p)2 

(CT  *  S)  +  l(CD  '  CTP) 

f-p£np (2-p) 

L  d-p)2 

(90) 


Now  (90)  is  larger  than  the  right  side  of  (89)  if 

CD  >  CTP  +2 (CT  -  s) 

Therefore,  if  CQ  is  large  enough  that  (91)  is  satisfied  and 
is  small  enough  that  (89)  is  satisfied  then  an  on-site 
system  is  optimal.  If  CD  does  not  satisfy  (91) ,  we  must 


(l-p) £np 
£np (2-p) 


(9 
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compute  the  right  side  of  (88)  for  several  values  of  N  to 
see  which  value  minimizes  it  and  then  compare  CT  with  the 
resulting  minimum.  If  CL  does  not  satisfy  (88)  or  (89) 
then  the  on-site  system  is  not  preferred  over  the  direct  de¬ 
livery  strategy  identified  as  Alternative  1  and  the  latter 
should  then  be  considered  for  use. 

The  initial  step  in  this  consideration  is  to  deter¬ 
mine  the  optimal  N  value  which  minimizes  equation  (15) . 

Next  we  must  consider  the  effect  of  the  constraints  of  time 
standards  and  truck  capacity. 

In  the  derivations  of  the  direct  delivery  models  the 
time  standard  T  was  assumed  to  be  measured  in  periods.  In 
reality,  it  is  measured  in  hours  or  days.  Therefore,  given 
a  time  standard  and  the  production  schedule  we  can  convert  the 
actual  time  standards  into  equivalent  periods. 

We  can  combine  the  time  constraint  given  by  (27)  and 
truck  capacity  constraint  as  follows: 

N  min  {M,  2T  +  1}  ,  (92 

where  M  is  the  truck  capacity  in  units  of  the  repair  part 
in  question. 

If  the  on-site  system  was  not  preferred  based  on  (88) 
or  C89) ,  then  the  optimal  N  determined  from  minimizing  (15) 
must  be  tested  against  (92) .  If  it  does  not  satisfy  (92) ,  we 
select  the  largest  integer  value  of  N  which  does.  We  must 
then  test  to  see  if  satisfies  (88)  when  this  N  value 

is  introduced  into  the  right  side  of  the  inequality.  If 
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now  satisfies  (88) ,  the  on-site  system  is  preferred  over  the 
constrained  direct  delivery  strategy.  Otherwise,  we  use 
Alternative  1  under  the  constrained  value  of  N. 

Alternative  2 — The  approach  in  developing  the  trade-off  model 
for  Alternative  2  is  similar  to  that  described  above  for 
Alternative  1. 

The  general  conditions  for  which  OECP(K)  <  ECP(N), 
where  equations  (85)  and  (35)  respectively  apply,  is  given  by 
(93)  . 


cL  <  <CT  -  SK)  j 


K,,  , n-K  L  cd(k~1) 

P  U-p)  +  - * - 


(93) 


When  K  =  1,  (93)  reduces  to  the  inequality  given  by 


(94) 


CL  "  (CT  '  s) 


w\ 


(94) 


The  right  side  of  (93)  when  K  =  2  is 


(CT  -  2S)  =  ff 

n=2 


n-ll 
1 


p2(l-p)n"2  +  CD 

T 


_  {r  v  n-1^2  M  ,n-2  D 

-  (C^  -  2s)  l  - p  (1-p)  +  — 

n=2 


-  <CT  -  2S>  ife 


V 

P  „  = 


n=2 


p(l-p)n_1  +  CD 


(Continued) 
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-  (CT  -  2s)  -  P  "  I%!~tnP  *  (1-P))J  +X 

•  <CT  -  2S»  +  4?f]  +  X  • 


The  conditions  under  which  (95)  is  greater  than  the  right  side 
of  (94)  are  expressed  by  (96)  and  correspond  to  K  =  1 
minimizing  the  right  side  of  (93) . 


CD 


2p_ 


(1-p) 


C_[-5.np  -  (1-p)]  -  s[ (1+p) (-fcnp)  -  2  (1-p) 


(96) 


If  CD  satisfies  (96)  and  satisfies  (94)  then 

the  on-site  system  is  optimal.  If  does  not  satisfy  (96) 

then  we  must  find  that  K  which  minimizes  the  right  side  of 
(93) ,  evaluate  the  corresponding  value  of  the  right  side,  and 
then  test  for  the  value  of  CL  .  If  CL  now  satisfies  (93) 
then  an  on-site  system  is  optimal.  Otherwise,  direct  delivery 
under  Alternative  2  should  be  examined  by  first  determining 
the  optimal  value  of  K  being  that  which  minimizes  (35)  and 
then  checking  this  K  against  the  constraints. 

The  time  standard  and  truck  capacity  constraints  under 
Alternative  2  can  be  combined  as 

K  <  min{M,  2pT  +  1}  (97) 

where  M  is  the  truck  capacity  in  units  of  the  repair  part 
in  question  and  T  is  the  time  standard  in  periods.  The 
feasible  value  of  K  is  the  largest  integer  value  satisfying 
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(97) .  If  it  is  at  least  as  large  as  that  K  minimizing  (35) 
then  the  problem  is  solved.  If  it  is  not  as  large  then  we 
must  proceed,  as  described  for  Alternative  1,  with  introducing 
the  constrained  K  value  first  into  (93)  and  testing  for  CT  . 
If  CL  now  satisfies  (93)  then  the  on-site  system  is  preferred. 
Otherwise,  direct  delivery  under  Alternative  2  with  the 
constrained  value  of  K  is  optimal. 

Alternative  3 — The  general  conditions  for  which  OECP(N)  <  ECP(N) 
when  equations  (87)  and  (55)  are  introduced  are  given  by  (98) , 


r 

00  . .  ,  k-1 

V 

U 

CT  +  CD(N-l)[i^  p  +  1 

-  s[l  +  (N— 1)  p ]. 

e  p.n~g _ r 

k*l  N  +  k  -  1 

A 

When  N  =  1,  (98) reduces  to 


CL  <  (CT  -  s) 


•^np 


(99) 


When  N  =  2,  the  right  side  of  (98)  is 


CT  +  CD  -  s(l+p) 


r  P(l-O) 


k-1 


=  [C_  +  Cn  -  s(l+p)]  —2 — =■  [-£np  -  (1-p)  ]  (100) 

T  D  (1-pr 


The  conditions  under  which  N  =  1  minimizes  the  right  side 
of  (98)  are  obtained  by  comparing  (100)  with  the  right  side 
of  (99).  The  result  is: 


(98 
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(101) 


c 


(1-p)  +  p£np 
-£np  -  (1-p) 


+ 


s 


2 

p  -  2£np  -  1 
-inp  -  (1-p) 


Therefore,  If  CQ  satisfies  (101)  and  satisfies  (99) 

then  an  on-site  system  is  optimal.  Otherwise,  direct  delivery 
under  Alternative  3  is  considered  further. 

The  time  standard  and  capacity  constraints  for 
Alternative  3  cannot  be  combined  because  equation  (76)  for 
the  expected  delay  per  unit  is  not  a  simple  function  of  N. 

As  a  consequence,  the  unconstrained  optimal  value  of  N  for 
Alternative  3  must  be  tested  against 


N  <  M 


and 


i  (N-2)  +  <  t  . 

/m  P 


If  the  constraints  are  satisfied  by  unconstrained  optimal  N 
then  the  solution  is  direct  delivery  under  Alternative  3  with 
this  value  of  N.  If  the  constraints  are  not  satisfied  then 
we  must  select  the  largest  integer  value  for  N  which  satis¬ 
fies  the  constraints  and  re-examine  (93) .  If  satisfies 

(93)  an  on-site  system  is  preferred.  Otherwise  we  use  direct 
delivery  with  constrained  N. 

The  Composite  Model — The  goal  of  the  trade-off  analyses  is 
the  determining  of  the  best  strategy  when  all  three  direct 
delivery  alternatives  are  considered  together  (the  composite 
problem) .  This  is  done  by  first  obtaining  the  best  solution 
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for  each  alternative's  trade-off  model  and  then  comparing 
these  solutions. 

Obviously,  if  an  on-site  system  is  optimum  under  ail 
three  alternatives,  then  an  on-site  system  is  the  solution 
to  the  composite  problem.  If,  however,  direct  delivery  is 
best  in  all  three  cases  then  the  optimal  solution  to  the  com¬ 
posite  problem  is  that  direct  delivery  alternative  which  pro¬ 
vides  the  lowest  optimal  total  expected  costs  for  the  given 
CT  ,  CD  ,  and  p  values. 

Now  suppose  that  one  alternatives  trade-off  model 
gives  an  on-site  system  as  optimal  but  another  does  not. 

When  the  two  results  are  compared,  that  alternative  favoring 
direct  delivery  is  automatically  preferred  since  it  was  pre¬ 
ferred  over  an  on-site  system  when  it  was  considered  by  itself. 
As  a  consequence,  if  only  one  direct  delivery  alternative  was 
found  optimum  under  the  individual  trade-off  models,  it  is  the 
optimal  solution  to  the  composite  problem.  If  two  trade-off 
models  favor  direct  delivery  and  the  third  does  not,  then  it 
follows  that  the  optimal  solution  to  the  composite  problem  is 
obtained  from  comparing  the  minimum  expected  total  costs  per 
period  for  the  two  models  favoring  direct  delivery. 
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VII.  SUMMARY  AND  RECOMMENDATIONS 


Summary — Trade-off  models  for  deciding  on  where  to  place  in¬ 
ventories  for  a  given  repair  part  needed  by  a  NARF  for  a  cer¬ 
tain  aircraft  or  component  rework  have  been  developed  in  the 
preceding  chapters.  Two  locations  for  the  inventory  were 
considered;  at  the  NSC  and  at  the  NARF.  Splitting  of  the 
inventory  between  the  two  locations  was  not  considered  in 
keeping  with  the  current  philosophy  of  the  NSCs  towards  NARF 
support.  When  the  location  was  assumed  to  be  at  an  NSC, 
three  direct  delivery  alternatives  were  considered  and  included 
both  scheduled  and  unscheduled  delivery. 

Expressions  for  the  total  expected  costs  per  period 
were  derived  for  all  alternatives  to  provide  a  basis  for 
comparison.  The  total  cost  was  assumed  to  consist  of  a  delivery 
cost  and  a  production  delay  cost  for  direct  delivery.  For  the 
on-site  inventory  system  at  the  NARF,  the  costs  were  assumed 
to  consist  of  a  labor  charge  and  a  delay  cost.  Constraints 
were  also  developed  to  reflect  the  impact  of  time  standards 
and  delivery  truck  capacity. 

Unfortunately,  the  complexity  of  the  various  expected 
cost  expressions  did  not  allow  for  optimal  solutions  to  be 
derived  analytically.  Therefore,  algorithms  were  developed 
for  using  the  models  to  resolve  the  question  of  where  to  site 
the  inventory . 

Recommendations — Further  understanding  of  the  models  should 


be  obtained  through  parametric  studies.  Some  of  these  are 


well  underway  and  suggest  that  between  the  direct  delivery 
alternatives  the  optimal  expected  costs  per  period  differ 
little  in  value.  Additional  studies  are  however  needed  and 
are  being  planned. 

As  more  of  an  understanding  of  the  models  is  obtained, 
refinements  will  undoubtedly  appear  appropriate.  In  fact, 
the  chapters  addressing  the  direct  delivery  models  contain 
some  preliminary  refinements  in  the  expected  delay  costs  ex¬ 
pressions  which  were  motivated  by  the  parametric  studies  which 
have  already  been  done.  Additional  issues  relative  to  delivery 
costs  for  scheduled  versus  unscheduled  delivery  have  also  been 
raised  in  Chapters  II  and  IV.  The  question  of  what  is  an 
appropriate  delay  cost  to  assume  for  an  on-site  system  also 
needs  further  study. 

As  was  mentioned  earlier,  the  models  have  been  restricted 
to  a  given  repair  part  for  a  certain  production  line  at  the 
NARF.  Expansion  of  these  models  needs  to  be  done  to  include 
multiple  sources  of  demand  within  the  NARF  for  the  parts. 

Finally,  an  interesting  additional  expansion  which 
also  seems  appropriate  is  the  case  where  two  or  more  repair 
parts  are  forced  to  have  the  same  decision  variable  value  in 
direct  delivery;  in  particular,  the  same  time  between  scheduled 
deliveries.  This  would  correspond  to  a  more  realistic  sched¬ 
uled  delivery  scheme  than  one  designed  for  each  part. 
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